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Abstract. We study the Cauchy Problem for the rclativistic Boltzmann equation with near 
Vacuum initial data. Unique global in time mild solutions are obtained uniformly in the speed of 
light parameter c > 1. We furthermore prove that solutions to the relativistic Boltzmann equation 
converge to solutions of the Newtonian Boltzmann equation in the limit as c — > oo on arbitrary time 
intervals [0, T] , with convergence rate for any e G (0, 2). This may be the first proof of unique 

global in time validity of the Newtonian limit for a Kinetic equation. 
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1. Introduction. The dynamics of a relativistic Gas is modelled by the rela- 
tivisitic Boltzmann Equation, which can be written as 

P^9^/ = C(/,/). 

The relativistic Boltzmann equation is the central equation in relativistic collisional 
Kinetic theory. Its collision operator is written in the physics literature [HITS] as 

C{Lh)^U ^[ ^[ ^W{p,q\p\q'mp'Mq')~f{p)h{q)]. (1.1) 

The kernel W{p, q\p', q') is called the transition rate. It takes the form 

W{p,q\p\q') = sa{g,9)S'^^+'\p^ + q'' - p^' - q^'), 

where a{g, 9) is the differential cross-section or scattering kernel; it measure's the 
interactions between particles. The rest of the notation is defined just below. 

1.1. History. First we will discuss a small part of the history of this fundamen- 
tally important model. Lichnerowicz and Marrot |46| are said to be the first to write 
down the full relativistic Boltzmann equation, including collisional effects, in 1940. 
In 1967, Bichteler [10] showed that the general relativistic Boltzmann equation has a 
local solution if the initial distribution function decays exponentially with the energy 
and if the differential cross-section is bounded. Dudyhski and Ekiel-Jezewska |20) , 
in 1988, showed that the linearized equation admits unique solutions in LP'. After- 
wards, Dudyhski |21| studied the long time and small-mean-free-path limits of these 
solutions. Then, in 1992, Dudyhski and Ekiel-Jezewska [22| proved global existence 
of large data DiPerna-Lions renormalized solutions |18| . In this paper, they use the 
causality of the relativistic Boltzmann equation, a result which they had previously 
established [23l[24]. This result has since been extended in [42ll43] . 

In 1991, Glassey and Strauss [28] studied the collision map that carries the 
pre-coUisional momentum of a pair of colliding particles into their momentum post- 
collision. In 1993, Glassey and Strauss [29] proved existence and uniqueness of smooth 
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solutions which are initially close to a rclativistic Maxwellian and in a periodic box. 
They also established exponential convergence to Maxwellian. In 1995, they extended 
these results to the whole space case [3^ where the convergence rate is polynomial. 
Analogous results with reduced restrictions on the cross section were more recently 
explored in [38] using the energy method from [32H35] . Very recently, the global sta- 
bility for full range of soft-potentials (in the sense of the assumption in [20]) and also 
rapid convergence in a periodic box has been established by the author in [59j . 

In 1996, Andreasson [T] showed that the gain term is regularizing and the strong 
compactness. This is a generalization of Lions [TTIHS] result in the non-relativistic 
case. In 1997, Wennberg [62] proved the regularity of the gain term for both the 
relativistic and non-relativistic case in a unified framework. In 2004, Andreasson, 
Calogero and lUner [5] showed that removal of the loss term for the Boltzmann equa- 
tion (relativistic or not) leads to finite time blow up of a solution. 

In the same year Calogero [TT] proved existence of local-in-time solutions inde- 
pendent of the speed of light and established a rigorous Newtonian limit on that time 
interval. Glassey, in 2006, established the first result proving global in time existence 
of unique mild solutions to the relativistic Boltzmann equation with initial data near 
the Vacuum state [27] for certain cross sections. Earlier near Vacuum results for the 
Newtonian Boltzmann equation were already obtained in 1984 by Illner and Shinbrot 
[39] using in part the iteration method from 1978 by Kaniel and Shinbrot [IJ; for 
further Newtonian results near Vacuum see also [5l[6j[T3l[37l[53l[60] and the references 
therein. Conditional asymptotic completeness for the relativistic Boltzmann equation 
was shown in 2007 [55]; the existence theorem from [37] does not allow the kind of 
decay which would be needed for the result of [55] . It is these three results [TTl[?7l[55] 
on the relativistic Boltzmann equation that our theorems below can extend. 

1.2. Notation. Before discussing our main results, we define the problem pre- 
cisely. The momentum of a particle is denoted by p^, fi ~ 0,1, . . . , N ioT N > 2. (The 
most important physical case is of course N = 3.) Let the signature of the metric 

be ( h • • • +). We set the rest mass for each particle m = 1. The momentum for 

each particle is restricted to the mass shell p'^p^i — — c^, with p° > 0. The speed of 
light is a universal constant c = 299 792 458 m/s in free space. In this paper we 
study the Newtonian approximation to Einstein's theory of special relativity, and so 
we vary the speed of light as a parameter c > 1. Further with p G M.^ , we may write 
p'' — {—p^,p) and similarly = {—q'^,q). Thus the energy of a relativistic particle 
with momentum p is po = \/ (? -\- The Lorenz inner product is then 

= -Pogo + P • (J- 

We write the standard inner product for vectors in asp ■ q^ 'l2iLiPi9i- 
We will now define a few important quantities: 

S -{p^ + g^')(PM + 9j = 2 {-p^q^ + c2) > 0. (1.2) 

Above s is the square of the energy in the centcr-of- momentum system: p + q = 0. 
The relative momentum is denoted 

9 = = ^J2{-P''% - c')- (1-3) 

Notice that s ~ + Ac? . This notation is used in [TS] , but it is different from some 
more modern authors by a constant factor. 
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The angle 0, which is the scattering angle in the center-of-monientum system is 
defined by 

cose = {p^-q^){p'^-cJ^)/g\ (1.4) 

For general four- vectors p^ , ,p^' , q^' the quantity on the r.h.s. could be unbounded; 
however when we assume the collisions are elastic 6 can be shown to be well defined 
|26[ p. 113]. Momentum conservation for elastic collisions is expressed as 



+ g'' = p^' + g^'. (1.5) 



One may also write (|1.5p as 



p + q = p' + q', 

where the first line represents the principle of conservation of energy and the second 
line represents the conservation of momentum after a binary collision. Notice that 
(|1.5I) and (|1.2I) together imply that s is a collision invariant: s{p^,q^) = s{p^' ,q^'). 
The Lorcntz inner product and g are also: p^q^ = P^'l'fj, and g{p^,q^) = g{p^' ^q^')- 
We see that the "transport term" is a Lorcntz inner product: 

^podt +p-Vx- 

The usual way mathematician's write the rclativistic Boltzmann equation is then 

dtf + p-^f=Q{fJ). (1.6) 

Here in comparison with the notation at the top of this paper, Q(/, /) = C(/, f)/po- 
This has become the contemporary standard notation in the mathematics literature. 

Above we consider / = f{t, x,p) to be a function of time t E [0, oo), space x £ M.^ 
and momentum p G M.^ . The normalized velocity of a particle is denoted 

P P n -7^ 

p — c — = — =. (1.7) 

Po ^/T+WJ^ 

Steady states of this model are the well known Jiittner solution's, also known as the 
relativistic Maxwellian. We write the normalized rclativistic Maxwellian as 

w X do. exp i-cpo) 

" 2(2^)(^-i)/2cif2(c2)" ^ ' 

Here if 2 is a bessel function, see [49l p. 449]. Standard references in relativistic Kinetic 
theory include [HIIllEaEelET] . 

Before discussing our main results, we will explain two different ways to reduce 
the collision integrals in ()1.1|) . These are analagous to the two well known expressions 
for the post-coUisional velocities in the Newtonian Boltzmann theory. Each expression 
that we discuss below converges in the Newtonian limit to the corresponding Boltz- 
mann equation in either the so-called uj or a notation from |63| . The first expression 
below has become the coordinates that are used in the vast majority of mathemat- 
ically oriented papers in the field. The second set of coordinates is not new to the 
physics literature, but it has not been widely used mathematically. These are the 
coordinates we use to prove Theorem l2.11 they make some of our estimates tractable. 
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Expression from Glassey-Strauss (1993). Glasscy and Strauss illustrated in 
[29| that a reduction of the collision integrals can be performed, without using Lorentz 
Transformations, to obtain 

Q{f,h)^ f lC,{p,q,Uj)[f{p')h{q')- f{p)h{q)]dujdq, (1.9) 

where the kernel is ICc{p,q,i^) = ^^^j^B{p,q,uj) with 

(po + qofpaqa 

B{p,q,uj) = c- 



[(p, + q,y-(^u;-[p + q]y]' ' 
In this expression, the post-coUisional momentum's are given by 



p' =p + aip,q,^)^ 

q' ^ q - a{p,q,u:)uj, (1.10) 



where 

KPo + qo)poqo{c^-{i;-^)} 



a{p,q,uj) 



1 2 



(po + qo)^ - • (p + q)} 
And the energies can be expressed as p'q = po + Nq and q'o = qa — Nq with 

dof 2uj- {p + q){pa{uj ■ q) - qo{uj ■ p)} 



No 



{po + qo? - {w • (p + q)y 



This last expression for the post-collisional energies can be derived from the mass-shell 
condition, see for instance [Ml p. 18]. These expressions clearly satisfy the collisional 
conservations (jl.Sp . The angle is then defined by plugging these into (|1.4p . The 
Jacobian [55] for the pre-post collisional change of variable (p, g) — >■ {p' ,q') is found 
with a non-trivial calculation whose result is 

d{p',q') ^ Ppgo 
dijp.q) Poqo' 

We do not use these coordinates in the proof of Theorem 12.11 in particular because of 
A^o above. Simply put there are too many angles, cj, in the expression for A^o, which 
makes it more difficult to utilize these expressions in a few of our proofs below. 

Now these variables are analogous to non-relativistic Boltzmann variables, the 
so-called co representation: 

p' ~ p + uj ■ {q — p)uj 

q' = q-uj ■ {q-p)uj. (1.11) 
With these post-collisional velocities, the non-relativistic Boltzmann equation is 

5t/+p-V/= Q(/,/), (1.12) 
with the collision operator given by 

Q{f,h) = / Koo{p,q,o^)[f{p')h{ii) ~ f{p)h{q)]duodq. 
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The kernel above is 



/Coo(p,q,w) \u- {p- q)\aoo{\p- q\,0). 



(1.13) 



Indeed when CTc Coo the formal Newtonian limit, c oo, of p.9p with variables 
(jl.lOp is this standard Newtonian Boltzmann equation. 

There is another expression for the collision operator, which we now elaborate. 

Center-of-Momentum Collision Operator. Another method of reduction is 
described in the physics literature |15j . which uses Lorentz Transformations in the 
center- of-momentum system to reduce the delta functions and obtain 



Qif,h) 



JV„HN-1 



Vc <Jc{g, e) [f{p')h{q') - f{p)h{q)] dwdq. (1.14) 



where Vc = Vcij), q) is the M0ller velocity given by 

Vc = Vcip,q) 




p_ 


_ q_ 


2 
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4po9o 

The post-collisional momentum in the expression (jl.l4p can be written: 

p + q 



P = 
I 

q = 



2 

p + q 



9 f , / iv , Ap + q)-^ 



inhere 7 = {pq + qo)/\/s. The energies are then 



(1.15) 



(1.16) 



2^ 
9 



LU ■ (p + q) 
w ■ (p + q). 



These clearly satisfy (|1.5p . The angle further satisfies cos 6 = k ■ ui, where fc is a 
unit vector. In fact, (|1.16p can be defined more generally only up to a (non-unique) 
Lorentz Transformation, A, satisfying (|5.16p as 



p" 
P' 



-A- 



q 

q' 



-A- 



—guj 



(1.17) 



This will work for any of the Lorentz Transformations in the appendix to thispaper; 
a complete mathematical derivation of these coordinates can be found in |58jrl 

These arc the rclativistic analogoue of the so-called <j representation for the New- 
tonian Boltzmann collision operator. In the non-rclativistic case, these variables are 



P 



P 



1. 



q\uj 



q = 



^Onc sometimes find's in the contemporary mathematics literature the collision operator l|1.14|l 
being used in conjunction with the previous coordinates l ll.lOI I. When things are done this way the 
cross section, a{g,8), would be implicitly redefined to be different from the cross section originating 
in the transition rate of the expression for the collision operator HI. 111 . 
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And the collision operator interaction is 

Qif,h) = i / \p-q\ a^{\p-q\,e) [f{p')h{q') - f{p)h{q)]diodq. 

Indeed, when (Tc a^o the formal Newtonian limit, c — > cx), of (|1.14p with variables 
(|1.16p is this system. One of the main results of our paper is to show that this 
Newtonian limit can be rigorously established on any time interval. 

2. Statement of the Main Results. We consider the relativistic Boltzmann 
equation (|1.6p in its mild form 

f*it,x,p)^foA^,p)+ f ds Q*{f,J,)is,x,p). (2.1) 



We are using standard notation to express a solution along characteristics succinctly 
as ff{t,x,p) = fc{t,x + pt,p). The relativistic collision operator Qdfc, fc){s,x,p) 
represents either (|1.9p or (|1.14p . and it will be useful to have the index c. 
In the Newtonian theory, the mild form is given by 

f*it,x,p) = foix,p)+ f ds Q*ifJ)is,x,p). (2.2) 



Here the characteristics have a slightly different formula: /^(t,.T,p) = f{t,x +pt,p). 
Additionally, the collision operator Q{f, /)(s, x,p) is here represented as in (|1.12p . To 
motivate the following relativistic developments, we will now review the usual solution 
spaces for a near Vacuum solution to (|2.2p . 

Solutions to the Newtonian equation are measured in the following norm 

ll/*l|oo = SUp^^^(^. 

t.x.p Poo [X^ p) 

Here the weight that we use is standard 

Poo{x,p)=exp{-a\x\^) ^l^ip), a,/3>0. (2.3) 

The non-relativistic Maxwellian is given by 

fiiv) = (27r)-^/2e-l"l'/2. 

Resuhs such as I5l|6l[13l|3ai39llll[53l[60| show that if < faix,p) € C°(Rf x M^) 
initially and furthermore there exists a positive number such that if 

< bo, (2.4) 



Pooix,p) 



then a unique non-negative global in time solution f{t,x,p) to the mild form of the 
Cauchy problem (|2.2p exists, and satisfies for all time the bound H/'^^Hoo < bi for 
a positive constant bi. We provide analogs of these statements for the relativistic 
Boltzmann equation in the following developments. 

In order to use the dispersion of solutions along the relativistic characteristics of 
the transport operator, x + pt, and thereby control the time integral we introduce the 
following weight, which will scale in the appropriate way as c becomes large: 

Pcix,p) ^ exp (-apokl Vc) J^P), a, /? > 0. (2.5) 
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Above J^{p) is the relativistic Maxwellian (|1.8p to the power /3. This weight function 
is motivated by the following invariant quantity: 

q'o C P'o C 

= c^t + ^|^|2 + + ^ |x + i (j5 - q)\' . (2.6) 
Po c qo c 

We found this invariant, or re-found it to discover, during discussions with several 
experts, that it has existed in the folklore for some time. It was furthermore re- 
cently recorded in Jiang |41j . This invariance follows from the first order energy and 
momentum (jl.Sp . It follows by expanding all the terms in the following way 

^ 1^ + tp\^ ^ P^|a;|2 _ 2tx ■ p + Cpot' - -c\ 
C C Po 

We will explain more towards the end of this section; the use of this invariant (|2.6p 
requires us to introduce a cut-off (|2.9[) in the cross section. The identity (|2.6p forms 
one of the crucial elements of our global existence proof of Theorem 12.11 below. 
We will use the following norm to measure the size of solutions to (|2.ip : 

!l/fl|c=SUp^^. 

t,x,p Pc{x,p) 

The supremum above is over t > 0, x,p G S.^ . Throughout the paper we will use 
the letter A to denote a generic positive constant, which is independent of the speed 
of light c > 1. These constants, A, may change from line to line and they will never 
depend on important parameters in the problem at hand. 

Hypothesis on the collision kernel: We assume the collision kernel in (|1.14p 
satisfies the growth/decay estimates 

0<(7e(c^,P,<z) < 1^1 (^1+ (^YT^) ^ +A2 .g-^l a(L.). (2.7) 

Above Ai and A2 are non-negative constants and o'(w) < f 1 < 00. We allow ai > 0, 
< 7 < —3. This includes the hard ball assumption (|5.19p as a special case. We also 
restrict the collision kernel to be supported on (|2.9p . 



See jl9| for a physical discussion of general assumptions, which include our con- 
ditions. It would also be possible to prove our theorems with some mild singularities 
in the angle o'(w). We are ready to state our main theorems. 

Theorem 2.1 (Uniform Existence). Consider the mild form of the relativistic 
Boltzmann equation ()2.ip with cross section of type ()2.7p supported on p.9p . Choose 



initial values < fo,c{,x,p) G C°(R^ x M^) such that 



p 

Pc{x,p) 



< b. 



There is a positive number bo with the property that if b < bo then a unique non- 
negative global solution fc{t,x,p) to the mild form of the Cauchy problem exists. 
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This solution satisfies the global in time estimates 



Pc{x,p) 



<\\f*\\c<b. 



The norms and of course the solution here are highly dependent on the speed of light, 
but the constant 's 60 and hi are explicit and do not depend upon c > 1 . 

Theorem 12.11 provides a proof of global existence for unique mild solutions to the 
relativistic Boltzmann equation (|2.ip for any value of the speed of light c > 1. Next 
we will study the Newtonian Limit of general solutions to (|2.ip . This next result 
applies to the solutions in Theorem 12.11 but it also applies under much more general 
conditions. We consider the space 

dcf 



II/cIIlil- = / dp ||./c(p)||l=o(ir«). 

It is in this norm that we measure convergence to the Newtonian Boltzmann equation 
()1.12|) in mild form ^2.2\i . For h G M.^ define the translation operator as usual by 
Tfif{x,p) = f{x + h,p) and similarly define by T,^/(x,p) = f{x,p + h). We have 

Theorem 2.2 (Newtonian Limit). Choose the initial data /o,{/o,c} € L^L"^ 
for all c > I which initially converge for some Ai > as 

II./o,c-/o||lilsc <^l/c^ 3fce (0,2]. 

Suppose these initial data {/o,c} each lead to global in time unique mild solutions 
to p.l|) . denoted fc{t,x,p), which satisfy the uniform estimate \\f^\\c < ^2 < 00. 
And similarly, suppose fo leads to a global in time mild solution to (j2.2p satisfying 
11/"^ 1 1 00 < 62 < 00. Here 62 need not be small. 

Consider collision kernels, {cTc} and (Too, which satisfy \crc\ ^ ctq < 00, and 

IfTc-'Tool < A2(l + |p|" + |(7r")/c^ 3m,n>0, fcG(0,2]. 

Further suppose /Coo(p, 9, i^) from (|1.13p is Lipschitz continuous in p, which is a con- 
dition on (Too • fule out concentration by assuming that for the limit 

!|r,?/o - fohlLr + \Kjo - fohlLr < Mhl \h\ < 1. 

Then for any fixed T > (which is allowed to be large) and < t < T , the solutions 
fc to (|2.1|) and f to (j2.2p corresponding to these initial data converge as: 

\\fcit)-fit)\\LiL^^ <A{5,T)/c''-', fce(0,2], V6e{0,k). 

All the positive constants aQ,b2, Ai, A2, A3, A{S,T) are uniform in the speed of light 
parameter c as needed, but all of these constants are allowed to be large. 

Remark 1. Of course, in Theorem \2.S\ if initially the {/o,c} satisfy the assump- 
tions of Theorem \2.1\ then the initial data convergence condition just stated guarantees 
that (|2.4p is satisfied for fo and the global existence of the Newtonian system (j2.2p is 
known, as long as the constant in (|2.4p is sufficiently small. 
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Additionally, the assumption \ac\ < <7q < oo is not crucial to Theorem \2.'A In 
fact singularities such as those in (j2.7|) could be allowed by modifications to proofs in 
this paper. However we leave these issues for future investigations. 

Remark 2. Previous works on Newtonian Limits for unique mild or clasical 
solutions to Kinetic equations, e.g. [3ll4l llll[T2lll6[|45[l54[|55j . are local-in-time in the 
following sense. They generally prove that there exists a time interval [0,T] upon 
which the solution exists independent of c. Their convergence proofs to the Newtonian 
approxim,ation in general do continue to hold validity as long as the solutions may 
exist. Yet these results do not provide a uniform estimate which establishes global 
in time existence of a unique solution as the speed of light becomes arbitrarily large. 
We show with the proof of Theorem \2.1\ that such estimates are available for the 
rclativistic Boltzmann equation. Then in Theorem \2. 2\ we establish that the solutions 
form a Cauchy sequence on any time interval. 

There are also Newtonian Limit results in the context of weak solutions [TliH] ; in 
particular ^ holds in the time periodic case and [7] establishes weak convergence to 
the Newtonian approximation globally in time for the Vlasov-Darwin system. 

Remark 3. Previous results on the Newtonian limit for the rclativistic Boltz- 
mann equation were obtained by Calogero in 2004- In this result, which holds 
for the periodic case , a local existence theorem is proven on an interval which is 
independent of the speed of light and then convergence is shown in L^L'^ upon the 
time interval in which solutions are known to exist. The proof uses the existence of 
solutions to the limit system, the Newtonian Boltzmann equation with Hard-Sphere 
kernel, which is known since |39[|44| . No convergence rate is given. 

Remark 4. A priori asymptotic completeness results in L^ are shown for the 
rclativistic Boltzmann equation by Ha, Kim, Lee and Noh in |36| conditional upon 
global existence of solutions in a functional space with uniform space-velocity decay. 
In the notation of this paper, one of the main results in |36j shows that given a solution 
fc{t,x,p) to ()2.1|) there exists a unique scattering state f^^c such that 

lim ||/#(i)-/+,cWllLi,=0. 

Theorem \2.1\ provides the first existence theorem with the needed uniform decay. 

Remark 5. In 2006, Glassey |27j gave the first proof of global existence for 
mild solutions to the rclativistic Boltzmann equation (j2.ip near Vaccum. This paper 
appeared 22 years after the corresponding result |39| for the Newtonian Boltzmann 
equation first appeared. Let us point out that a crucial step in |39| is the following 
second order Newtonian symmetry for the variables (jl.lip ; 

\x + tp\^ + \x + tq\^ = \x + tp'\^ + \x + t(if , 
which follows from the Newtonian conservation of energy and momentum: 

\pf + \qf = \p'f + \q'f 
p + q=p' + q'. 

It was noted in |27| and 10 years earlier in |26j . that this symmetry fails in the 
rclativistic case; which is true because in particular now the energy and momentum 
are both first order (jl.Sp . Glassey 's proof proceeded otherwise by introducing a new 
weight function 



Pg{x,p) = e-"P«(l + |a; X p|2)-(i+'5), < 5 < 1. 
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Rather than having a cancellation of the post-collisional momentum via symmetry, 
the weight function pg cancel's very nicely the worst effect of the post-collisional mo- 
mentum because of orthogonality. The condition 5 < 1 is used to generate the usual 
reproducing property of the non-linear estimates for the weight function. 
The set of cross-sections allowed by this result are given by 

\uj ■ [qx p)\d{uj) 



d^^-±- < c\z\-\ 

This assumption in particular does not include the hard-ball case (j5.19p . 

On the other hand, Glassey's result |27j does not require the assumption which 
we explain just now. In the proof of Theorem 12.11 we cut-off part of the angular 
integration domain S^^^ as in (|2.9p below. To describe this we recall the identity 
(|2.6p . which is equivalent to the following 

— \x + tp\ H \x + tq\ = — Ix + tg ^ \x -V tp \ 

c c c c 

P'n Pa 10 

Let us carry out this computation more completely, making use of (jl.Sp . 

P^\^ + tpf = P^\x\'-2tx-p+^\p\H' 
c c c 

= — \x\'^-2tx-p + cpot'^- —c^ 
c Po 

= ^|x|2-2te-p' + cp't2--c3 
c Po 



c3 



c c p'a Po 

c Po Po 

This shows how to obtain the invariance (|2.6p . 

Since the terms of the order t^ in the parenthesis of the previous display do not 
have a sign, they are hard to control pointwise. Alternatively, the weak dispersion 
from the transport operator is unknown and unlikely to induce the kind of decay 
which would allow the use of this weight without the cut-off that we are going to 
introduce now. 

For a given constant B > and a number < a < 1 and t > 0. we define 

,/ N -B CtOn \x + 1 (p — q)\^ / c „ ,„ „, 

h,=.h{x,p,q,t,c) = -+a-^ ^ >0. (2.8) 

We remark that he can generally be quite large for various values of it's arguments. 
Now we define the cut-off set 
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As already noted, we require that the differential cross sections CTc are all supported 
on this set (|2.9p . We also point out that 

B a\x + t{p-q)f ^ 
he rioo = ^ + a > 0, as c — >■ oo. 

Furthermore 

Lemma 2.3. limc^oo^c = 

We will prove Lemma 12.31 in Section [3] Lemma 12.31 shows that the cut-ofF does 
not pose a restriction in the Newtonian Limit. 

This assumption may be considered as a relative of Grad's angular cutoff for the 
relativistic Boltzmann equation. However, in general, for some values of p and q this 
assumption will cut-ofF a larger proportion of S^~^ than the assumption of Grad. We 
would like to see a future theory which allows collision kernels such as (|2.7p and does 
not require this cut-ofF. In our opinion such a future theory may require either (1) the 
discovery of a new invariant, or (2) a substantially different method of proof from the 
original Illner-Shinbrot [31] approach. Our main motivation in this work to illustrate 
that one can prove a global in time Newtonian limit for this Kinetic equation. 

The rest of our paper is organized as follows. In Section [31 we prove a series of 
pointwise estimates on algebraic quantities which will be necessary in the following 
sections. In Section[4l we prove the global existence Theorem l2 . ll using the coordinates 
pTTej) . And in Section [3 we prove the global Newtonian limit Theorem 12.21 using 
instead the Glassey-Strauss coordinates (|1.10p . 

Then in Appendix A, we give examples of Lorentz Transformations that lead to 
the coordinate system (|1.16p which we use in Section [4| of this paper, as described 
earlier in ()1.14p . And finally in Appendix B, we describe some collisional cross sections, 
<T{g, 0), that can be found in the physics literature. 

3. Algebraic Estimates independent of c. In this section, we prove a series 
of pointwise asymptotic estimates that will be useful in the following sections. We 
first prove an estimate which is stronger than Lemma [2.31 In particular, this lemma 
allows us to show that the solutions from Theorem 12.11 satisfy the Newtonian Limit 
Theorem 12.21 and in the limit condition (|2.9I) disappears. In fact more is true 

Lemma 3.1. For a fixed collection p, q, T, B , with < t < T , there is a constant 
— c* (p, q, T, B) such that for c > c^ we have 

where lB^(a;) is the indicator function of the set Be defined in (j2.9p . 

Notice that Lemma implies Lemma [273] More generally, this estimate implies 
that the solutions from Theorem 12.11 satisfy the results of the Newtonian Limit The- 
orem 12.21 Let us be more specific. Consider a cross section (Tc which satisfies the 
conditions in Theorem 12.21 Then define Cc — ^c^b^- Now we may construct global 
solutions with this cross section, CTc, using Theorem l2.1l Further the final lower bound 
in the proof of Lemma |3. II will imply that Lemma 15.21 below remains valid for Cc and 
this is enough to prove the convergence property in Theorem 12.21 

Proof. We start with the condition on the complement of Be that 

1 1 1 1 < ^ 
Po Qo Po q'o ~ c3' 
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Multiply by Po'ZoPo'Zo obtain 

{qo + pa)poqo ~ Paqaiq'o + p'o) < -^poqop'oq'o- 

Thus 

, , ^ h PoqaPWo 

Poqo-Paqo< — o 7 ] r • 

c-^ (PQ + qo) 

We use the energies of the coordinates (|1.16|) 

, Po + go , 5 / , N 

Po = o + flc, Oc = TT-F^ ■ (P + q) 

Z Zy/S 

I Po + qo 

go = — 2 "'^^ 

Notice 

/ , f Po + qoY 2 Po , go , Pogo 2 
Poqo=[-^) -a^-- + - + — ^a^. 

We plug this in to obtain 

Po , go Pogo _ 2 ^ _ ^ PogoPogo 



This implies 



Let 



Equivalently 



4 4 2 c3(po + go) 



2 ^ ^ PogoPogo (PO - go)^ 



-a; < 



(po + go) 4 
|p + g| COS0 = w • (p + g). 



COS^e> 4 PogoPogp S ^ (Po-go)^ g dof ^2 

~ |p + gp c3 (po + go) 5^ |p + gp 5^ 

We now see that goes to infinity with rate O(c^), since s = + 4c^. 

We will derive a sequence of (sometimes crude) lower bounds for £^ which allow 
us to elaborate on this point. First of all there is still some weak angular dependence 
in £'^.. However we can bound from below 



2p;g; = cVl + |p'IVcVl + |g'IVc2 > cpo. 
We will prove this last inequality in Lemma [3T2] below. We therefore have 

2 h cpIqo s 



> 



\p + q\ (Po + go) 9^ ' 
For < t < T, h > Y- Since s = + 4c^, we may further bound 

2 B plqa 4 



t > 



|p + g| r (po + go) 
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Now 

,|2\ / l.,|2\l/2 



J \ (? ) V (? 



and 



Thus 



(po + go) 2 

Furthermore 17 < |p — g|. AU these estimates lead to the crude lower bound 

cos^ 9>ll> -—^ -^c\ 

This easily establishes the Lemma by showing Be = §^~^ for large c. □ 

We remark that this last lower bound is crude for large values of \p\ and \q\; in 
particular as can be seen directly from the quantity £^ above the set Be does not 
rule out large momentum in general. This Lemma 13.11 is our main crucial use of the 
coordinates (|1.16p . The rest of the proof of Theorem 12. II is coordinate independent. 
Now we give a series of pointwise asymptotic estimates. 
Lemma 3.2. We have the pointwise estimate 



p'oq'o = cVi + b'l VcVi + l?T/c2 > CPO- 

And further 

\p-q\ 



\p-q\> 9>c 



y/poqo 



Proof. The estimates in this lemma are minor refinements of estimates in |29j 
(which are done with c = 1). Our estimates allow control of the scaling in c > 1. 
By energy conservation (jl.Sp we have 



^1 + |p|2/c2 < ^l+|p'|2/c2 + ^l+|g'|2/c2. 

Squaring both sides 



1 + bl Vc' < 1 + bf /c^ + 1 + \q'\''/c' + 2Vl + b'IVcVl + l'7'IVc' 

<2(l + b'|Vc2 + l + k'|Vc') 

<4(l + b'|Vc^) {l + lqf/c') 

Taking square roots and multiplying by yields the claimed inequality. 
Next, by a difference of squares, we can write (|1.3p as 



V Poqo + p- q + c^ J V Poqo + p- q + c^ 

c^\p-q? + b X 



= 2 



Poqo + p- q + c^ 
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And the estimate follows easily. The upper bound proof is the same as in [52]. □ 
In the next Lemma, we estimate the pointwise asymptotic convergence rates of a 

few quantities that will need to be controlled in the sequel. 

Lemma 3.3. For c> 1 we have, under the assumptions of Theorem \2.'A that 

1 + IpI" + lol'" 

\lC,{p,q,u:)~lC^{p,q,Lo)\<A ^ 1^' ^ 1'^' , fce(0,2], 3m, n > 10. 

c 

Furthermore 

\JCoo{p,q,u:)\< Aao{l + \p\ + \q\) . (3.1) 
We have the upper hound for ICc as 

|/C.(p,g,w)| <^ao(l + H)(l + |9p)5/2. (3.2) 
By symmetry, since ]Cc{p,q,uj) ~ ICciq,p,(^), we have 

\IC4p,q,uj)\<A<Jo{l + \q\){l + \pff/^. 

Notice that all of the constants A > above are independent of c> 1 . 

Proof. These basic pointwise estimates were shown for the most part in Calogero 
[TT] . In particular the first estimate for the difference of \ICc{p, q,uj) — K-ooiPi q,i^)\ 
follows from [TTl Lemma 1 part (b)] when combined with the estimate for the difference 
[cc — tJool in Theorem 12.21 and (|3.1|) . Note that JCc in this paper is defined differently 
from /Cc in [llj . Furthermore, p.ip follows trivially from the assumptions in Theorem 
12.21 (in particular that |(Too| < <^o < oo) and (|1.13p . Lastly the estimate p.2|) is proven 
as an intermediate step in the proof of [11] Lemma 1 part (c)]. □ 

We also have estimates for the difference between the post-collisional relativistic 
momentum (|1.10p and the Newtonian post-collisional momentum (jl.lip as 

\p'~p'\ + \q' -q'\<A{\p\ + \q\f/c^ (3.3) 

This is shown in [ill Lemma 1 part (a)]. 

Furthermore, for the difference of normalized velocities we have 

<A(i±^). (3.4) 



This estimate is proven with a first order taylor expansion of the denominator. The 
first terms cancel and the remainder terms are bounded as the upper bound in p.4|) . 

Lemma 3.4. The relativistic Maxwellian J{p) as in (jl.Sp is bounded and has 
exponential decay independent of the speed of light. In particular 

A-ie-^lfl' < J{p) < Ae^IPl, 

where A > does not depend on the speed of light c. 

Proof. In (|1.8p K2{x) for say a; > 1 is the Bessel function: 



2 /"OO 
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ds. 



We translate s s + 1 to obtain 

3 Jo 

And we change variables s — )- s/x to obtain 

1 /"Oo \ 3/2 

From the last representation, we conclude that (for a; > 1) 

A-^x-i/^e^^ < K2{x) < Ax-^^^e~\ 
Plugging this inequality into (jl.Sp we conclude 

A-^e''\-''P° < J{p) < Ae'^'e-^P". (3.5) 

We claim that 



X2(c) c2 - cVi + bl7c2 < X2(i) = 1 - yTTW, c>l. (3.7) 



xi(blVc') = 1 + - Vi + blVc^ > (3.6) 

2 



Then p.6p establishes the lower bound and p.7p establishes the upper bound. 
Wc first establish p.6p . consider Xi(^) for a; > 0. We differentiate 

Since Xi(a;) is increasing and xi(0) = we have p.6p . 

We finish off the proof by establishing (|3.7p . We will show that X2(c) < for any 
fixed \p\ > and for all c > 0. We compute 



Notice that X2(0) = ^|p| < and X2(c) is continuous. Suppose X2(c*) = for some 
c* > 0. Then 

2c*iic*r+\p\Y^'^{2ic*r+\p?). 

After squaring both sides we observe that this implies |p|^ = 0, which is a contradic- 
tion. Hence X2(c) < for any fixed \p\ > and for all c > 0. □ 
The following estimate (|3.8p will be important in Section [5j 
Lemma 3.5. As c f oo, we have the following pointwise convergence 

For any positive function h{c) < Awfc, there is a uniform constant A> Q such that 

J{h{c)) < Ae-'''(^)/2_ (3 8) 

By J(h(c)), we mean J{p) evaluated at \p\ = h{c). 
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The last part of this lemma p.8p is a subtle asymptotic estimate. Despite the 
fact that it is not hard to prove, this estimate is crucial to speed up the convergence 
in our arguments below for the new rapid convergence rate l/c?~^ in Theorem [ 

Proof. For < cc, we have the following Taylor expansion 

1 



VlT^ = I + -X + R{x) (3.9) 



Recall (|1.8p and note that for c large 

Using this, (|1.8p and p.9p we have 

J{p) - (2^)-^/2gC=-cpo ^ ^^^y 

This completes the proof of convergence to a Newtonian Maxwellian. 

We now prove the upper bound (|3.8p . We use the upper bound in p.Sp . First 
the taylor expansion p.9p yields 



2 9 9/ bp bp 9 /bp 



2 V c' 

Of course we are working on the set \p\ — h{c) < Ai^/c. Clearly on this region the 
expression for R in (|3.9p gives the bound 



R 



C2 



1 

< -. 



Plugging the last two estimates into the upper bound (|3.5p yields p.8p . □ 

Let r(c) = (logc)"^ for i < ai < 1. Given A > 0, e > there is a uniform 
constant B = B{e) > such that 

Furthermore there is a constant Bi{k) > such that 

^-Mr^ic)^^M^ Vfc>0. (3.11) 

These estimates are obvious by after taking the logarithm of both sides. 
In the next section we prove 

4. Global Existence uniformly for c > 1. In this section, we prove a uniform 
(in c > 1) non-linear estimate which is the key step in our uniform global existence 
Theorem 12.11 Global existence of unique mild solutions can be deduced from this 
estimate via a standard fixed point argument for small data. After the proof of 
Theorem 14.11 below . we outline how this works. 

We split the collision operator ()1.14p into its "Gain" and "Loss" terms as 



Q+(/, h)^ V, a{g, 9) f{p')h{q') dtudq 

S- (/, h)^ [ Vea{g,e) f{p)h{q) dtudq. 

JR«xS"-i 



Global Newtonian limit for the Rclativistic Boltzmann Equation 



17 



With this notation in hand, we prove the foUowing non-hncar estimate. 

Theorem 4.1. For collision kernels satisfying (|2.7I) and (j2.9p . we have 



ds 



Q*{f,h){s,x,p) <Ap4x,p) \\f*\\c\\h*\\ 



Similarly, for the loss term 



ds 



Qtif,h)is,x,p) <Ap,ix,p) \\f*\\c\\h*\\ 



The constants A > above are independent of the speed of light. 

This estimate may formaUy look very similar to the Newtonian estimates near 
Vacuum, e.g. [55]. However the main new features of this estimate are twofold. We 
have firstly made use of the invariant (|2.6p combined with (j2.9|) in order to exploit 
the weak dispersion along characteristics. Equally important, we have introduced the 
correct weight (|2.5p and shown that this allows us to obtain the estimate in a way 
which is invariant with respect to c > 1. Now we are ready for the proof. 

Proof. We begin with the loss term. By definition 



ds 



Q*{f,h){s,x,p) 



ds 



'xS«-i 



diodq Vc <ycf{x + sp,p)h{x + sp, q) 



ds 



JVxSN-l 



divdq Vc cTc f*{s,x,p) h*{s,x + s{p - q),q) 



< \\f*\\c\\h*\\c Pc{x,p) I ds I dq duj Vc <Tc{g,0) pc{x + s{p - q),q). 



We use p.5|) . Lemma 13^ and then p.7p to conclude 



ds 



dq 



ds 
t 



duj Vc cr(.g, 0) pc{x + s{p - q), q) 



dq 



duj Vc a{g,0) J'3(g)e-"9«l^+-^(P-9)l'/'= 



< 



ds 
dq 



dq / du! Vc a{g, 9) 
dw Vc <j{g,0) e-'^l9l 



-0\<l\„-O'\x+s{p-q)\'^ 



< 



dq 



ds e-"l^+^(P-«)l' 

duj a{g,0) e-'^l'?! < A. 



In getting to the last line above we utilized the standard estimate 



ds e-«l^+'*(P-9)l' < 



TT 1 



a\p-q\ 



This follows easily by completing the square in the exponent of the exponential; see 
for instance [55]. We have also used in the same spot the estimate 



Vc 



c 

\p-q\-2 



P_ _ ± 



1 



1 



\p-q\ 



18 



Robert M. Strain 



which fohows from (jl.lSp . In our estimates for the loss term we do not need to use 
the restriction into the set (|2.9p at all. 

Furthermore, we explain the last inequality which will follow from (|2.7p 



duj cr{g,0) e^'^l"! <A dq dio {l + g-'') e"'^!"! < A. 



dq 



The last bound follows easily from our independent of c bound on g in Lemma 13.2 
Next we consider the gain term, for which we have the estimate 



ds 



Q*if, h){s,x,p) 



= / ds 



dq duj Vc o f{s, x + sp, p')h{s, x + sp, q) 



ds 



<!l/*llc!|/^*l|, 



dq 

t 



duj VcU f*{s,x + s{p~p),p') h*{s, X + s{p~- q),q') 



ds 



dq duj Vc a pc{x + s{p - p'),p')pc{x + s{p - q'),q'). 

Next we use the weights pc and their arguments to control the remaining integrations. 
Using ([23]), dm) and dUS]) we have 

Pc{x + s{p - p'),p')pc{x + s{p - q'), q') = pc{x,p)pc{:x + s{p - q), q)e~^*' . 



1 

7o 



> -he 



On Be, in 1^3}, we have 

1 1 1 
— + T 

kPo qo Po 

Recalling (|2.8p . we have the bound 

Combining this estimate with the weight (|2.5p and recalling that < a < 1 we obtain 
the upper bound 

Pc{x + s{p - q),q)e-^'' < e^e-(i-'^)"*l"+*(P-«)l'/= J^(<;). 

We conclude 

Q*{f,g)is,x,p) 



ds 







< \\f*\\c\\h*\\c Pc{x,p) 



ds 



dqduj Vc <7c e-S-(i-'^)"9ok+t(p-g)lVc j/3(^) 



<A\\f 



*\\\\h*\ 



Pc{x,p). 



The final integral is controlled in the same way as for the loss term. □ 

Now that our non-linear estimate has been established, it becomes a standard 
exercise in well known and well exposited techniques, e.g. |39ll44j and for instance 
[26) . to prove global existence of unique positive mild solutions. In particular the fixed 
point argument and Kaniel-Shinbrot iteration scheme described in [27] can be used 
without any modification in combination with our scale invariant non-linear estimate 
in Theorem 14.11 in order to prove global in time existence of unique mild solutions 
with a small data condition which is uniform in the speed of light c > 1. 
In the next section we consider the Newtonian Limit. 



Global Newtonian limit for the Rclativistic Boltzmann Equation 



19 



5. Global Newtonian Limit as c oo. In this section, we work with the 
relativistic Boltzmann equation (|1.6p in the forward mild form 

fc{t,x,p) = Sc{t)fo,c{x,p) + ds Sc{t - s)Qc{s,x,p). (5.1) 

Jo 

We are using the definition Sc(t)f(x,p) = f{x — tp,p), where we recall from (|1.7p that 
p ~ c— ~ I ^ clearly depends upon c. Notice that this mild formulation can be 

easily interchanged with the backward formulation that we used earlier p.ip ; this one 
is however more convenient for the Newtonian limit. Existence of such solutions is 
covered by Theorem l2.1l Above we also use the abbrieviated notation Qc = Qc{fc, fc)- 
We also study the Newtonian Boltzmann equation (|1.12p in the forward mild form 

fit,x,p) = S{t)fo{x,p)+ f ds S{t-s)Q{s,x,p). (5.2) 

We are here also using the definition S{t)f{x,p) = f{x — tp,p). The existence theory 
for (|5.2p was explained in Section [21 For the limit system, we have continuity. 

Lemma 5.1. Let the conditions in Theorem \8.2\ hold. Given T > and any 
< S < 1 there is a constant A = A{S, T) > such that for any < t < T we have 

sup (ii</(t) - /(t)iiLiL~ + ii</(o - fmLiL^) < m', \h\ < 1. 

\u\<\h\ ^ " ^ 

Lemma 15.11 proves the continuity of the mild form of the Newtonian Boltzmann 
equation which will be an important part of our convergence theorem. 

Proof. To estimate the x difference we subtract the two equations (|5.2p as 

rUip) - .fip) = T^S{t)fo{x,p) - S{t)fo{x,p) 

+ / ds Tf^S{t ~ s)Q{s,x,p) - S{t - s)Q{s,x,p), 
Jo 

which can then be further expanded as 

= T^S{t)fo{x,p)~S{t)fo{x,p) (5.3) 

diodq /Coo {r^,S{t - s)f{p') - S{t - s)f{p')) Tf^S{t - s)f{q') (5.4) 



ds f dojdq /Coo S{t - s)f{p') {T-^S{t - s)f{q') - S{t - s)f{q')) 

ds f dujdqICoo {T^S{t-s)f{p)-S{t-s)f{p)) T^S{t^s)f{q) 

* ds f dujdq /Coo S{t - s)f{p) {T^S{t - s)f{q) - S{t - s)f{q)) . 

The LpL^ norm of (|5.3p is easily controlled 

MS{t)fo - S{t)fo\\LlLr = Mfo - fohlL^, 

which is simply a statement of translation invariance of the norm. 
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For the other terms, we spht into an unbounded and a bounded region 

1 = l|p|+l9l>fl + l|p|+kl<i?' (5-5) 

where we fix a suitable R > below. 

We take the supremum in x and integrate over p. On the unbounded region, 
\p\ + \q\ > R, we see that each of the terms ()5.4p and after are bounded above by 

62/ ds f doudqdp IC^ lipi+igi>u ^^^ip) ^^^iq) <b2aoTe-^^'/\ (5.6) 



We have used the uniform bound for solutions in Theorem 12.21 p.3p . and p.ip . 
On the bounded region we have the following upper bound for the term (|5.4p 

< f ds f dujdqdplipi+igi^nlCooMf{p')-f{p')\\Lr\hhfiq')\\L^ 

Jo JR«xR"xS"-i 

< ds dujdqdp lipi+igi<4ii ICoo\\tU(p) - f(p)\\L^\\'rhfi<l)\\L^ 

Jo JR"xR"xS"-i 

< ARao f ds f dwdqdp \\TfJ{p) - f{p)\\L^\\rl^f{q)\\L^ 

Jo jR«xR"xS"-i 

<ARj\s ||r;^/-/|U.ij.||/|U.io. 

<AR f ds \\T^f - f\\LlL^\\f*\\oo<Ah2R f ds \\TfJ-f\\LlL^. 

Jo Jo 
Above the constant A > is independent of R. We just used the pre-post collisional 
change of variables, p.ip and the assumption of Theorem 12.21 for the uniform bound 
in terms of 62- Notice that all the other terms after (|5.4p can be treated in exactly 
the same way. We have shown 

+ A2R f ds WJ-fhlL^. 

Jo 

By the Gronwall inequality 

\\rU - fhiLr < {\Kfo - fohiLr + AiTe-^«V4) ^A^^nr ^ 

Now choose R = (-log|/i|)"' for any ^ < ai < 1. Using (fO)]) and ([3lT|) with 
c = i/\h\ and \h\ < 1 we obtain for any small e > that 



WrLf fhlLr < A [\h\-lrUo foh^^Lr + I 

The proof of the estimate for ||t^/ — is completed by using Lipschitz property 

on the initial data in Theorem 12.21 

For the estimate of |1t^/ — /Wl'^l^, we write out the difference of (|5.2p as 

rlfip) - f{p) = rlS{t)fo{x,p) - S{t)fo{x,p) 

+ / ds Tf^S{t - s)Q{s,x,p) - S{t ~ s)Q{s,x,p), 
Jo 

= Ti + T2. 
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We can expand the first term as 



Ti = fo{x - t{p + h),p + h)- fo{x - tp,p) 

= {.fo{x - tip + h),p + h) - fo{x - t{p + h),p)} 

+ {fo{x-t{p + h),p)- foix - tp,p)} . 



Thus, 



lirilliiLoo < Wrfja - /o||lil~ + hthfo - MIl^l^ 

< - M\LlL^+T\\THf, - /oIIlil-. 

Here we used the translation invariance of tlic norm to pull the fixed constant out of 
the supremum. More precisely, for any integer n > 1 we used the estimate 



1=0 

n 

< \\^hfo - /oIIlil- < n\\T^fo - /oIIlil-, 

i=0 

which follows because of the translation invariance of these norms. Technically, we 
should have supposed that T was an integer, but this does not cause difficulties. 

The other term T2 involves the collision operator. In this spot we use the New- 
tonian variables (|l.lll) . We further split as 12=72 + ^3 with given by 

T2 = / ds / dujdq {TllCoo){p,q,uj) f{s,Xh,p'h) fis,Xh,q'h) 







t 

ds I dujdq ICocip,q,uj) f{s,x,p') f{s,x,q'). 







The remaining term is then given by 







73 = -/ ds dujdq (rllCac.) {p,q,i^) f{s,Xh,p + h) f{s,Xh,q) 







+ ds diLidq ICoo{p,q,i^) f{s,x,p) f{s,x,q). 



We have used the following simplifying notation: x^ ^ x — (t — s){p + h) and similarly 
X ^ X — (t — s)p. Also with (|l.lip we define 

Ph ^ P + h + uj ■ [q - p - h) uj, q,^^q-uj- {q-p-h)uj. 

For both ^^'^ 2^3j use the splitting (|5.5p . On the large region l|p|+|q|>/j, we 
again have (j5.6|) for both T2 and T3 (the small < 1 causes no disruption). 
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On the bounded region, l|p|+|g|<fl, we expand as 

^ I I ^^^^ ■'■|PI+I9I:^-R^°° {f{s,Xh,Ph) ~ f{s,Xh,p)} f{s,Xh,qh) 
+ J ds J dujdq l\p\+\q\<R K,^ f{s,Xh,p'){f{s,Xh,q'h.) ~ fiS:Xh.,q')} 
+ [ ds [ dujdq lipi+iyi<R ICac f{s,Xh,p){f{s,Xh,q) - f{s,x,q)} 

ds / dx^dq l|p|+|g|<_R /Coo f{s,x,q') {f{s,Xh,p') - f{s,x,p')} . 



As a result of the hypothesis in Theorem 12. 2[ since /Coo is Lipschitz continuous we 
have 

IK/Coo) - /Cool e-'^d^'l'+l'?!')/^ < A|/i|e-^(l^'l'+l«l')/8. 

Here there could be a smooth polynomial appearing in \p\ and/or \q\ which is controlled 
by the exponential decay. With this estimate, the first term in integrated over dp, 
has the following upper bound estimate 

J^s J dujdqdp i(tp/Coo)-/Coo| umh^ 

<AblT J diudqdp |(t,^/C) -/C| e-^^IPl'+l^l'^/^ < AblT\h\. 



The estimate above follows from the uniform bound for solutions in Theorem 12.21 and 
(|2.3p as well as the pre-post invariance of the exponentials. Notice that the translation 
of a small h does not cause difficulties. 

For the next terms in with the bounded region, we will use several times that 

\q',,-q'\ + \p',-p'\<m 

This is a trivial consquence of Now the estimates for the next terms in 

work as in a previous case. In particular 

t 





ds / dwdgdp l|p| + |,|<^ |/Coo| \\f{Ph)-fip')\\L^\\f{<lh)\\L- 



<A'fds f dujdqdp lipi+i,i^n \ICoo\ sup ||(tP/)(p')-/(p')IU.-//'(9') 

Jo J l«l<4|'i| 

<A' f ds /rfc^dgdpl|p|+|,|<4i? I/Cool sup ||(rP/)(p)-/(p)||L^//2(^) 

Jo J l"l<4|/i| 

<A'R f ds sup \Kf-f\\L^L^ 
JO l"l<4|/i| 

<AA'R f ds sup WrU^fWLlL^ 

Jo \u\<\h\ 
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Here we have used the estimate for solutions in Thcorcm l2.2l and p.ip . All of the other 
remaining estimates in are exactly the same as this one, or a previous estimate in 
the X variables given earlier in this proof. The loss term estimates arc the same but 
easier because they do not require a pre-post change of variable. 
For the T3 term we employ the following splitting 

- Ts = y ds J dujdq { (r^/Coo ) - A^oo } f{s,Xh,p + h)f{s,Xh,q) 

+ y ds J dcodq JCoo {f{s,Xh,p + h) - f{s,Xh,p)} fis,Xh,q) 
+ J ds J dujdq JCoo {f{s,Xh,q) ~ f{s,x,q)} f{s,Xh,p) 

ds J dudqKryo {f{s,Xh,p) - fis,x,p)} f{s,x,q). 

All of these terms have a corresponding expression in the expansion above (except 
there are no primed variables here, which simplifies things). Thus the estimates for 
T3 are exactly the same as the ones for Tj, and they are proven in the same way. 
Collecting all of these estimates, we have shown 



\rU - fW^Lr < Kh - MIlil^ + A'T sup IK/o - Ml^Lr + ^iTe 

\u\<\h\ 



+ AT\h\ + A^i? /* ds sup UtH - /lU.^oo + IK/ - IWlil^) . 

We apply the Lipschitz condition on the initial data in Theorem 12.21 as well as the 
prior estimate for Hr,^/ — to obtain for any small e > that 



\u\<\h\ P X V J 

+ A'R [ ds sup \\rPf~f\\L.r^. 
Jo \y-\<\h\ 

By the Gronwall inequality 

\u\<\h\ P X V J 

As before choose i? = (— log for any ^ < ai < 1. Using p.lOp and p. lip with 

c = l/\h\ and < 1 we obtain for any small 5 > the statement of Lemma [5. II □ 

We now give our main estimate for the difference of the two collision operators. 

Lemma 5.2. Assume the conditions of Theorem \2.2\ hold. Suppose that c > 1, 
i < ai < 1, and < t < T . Then for any e € (0, 2] we have 

mt)Q4t) - smit)u.L^ <Ai\ogcr ii/c - fu^L^ + 

These collision operators are defined as in (jl.9|) and (|1.12p . Furthermore, as usual, 
the positive constants A, A{T) are uniform in c. 
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Proof. From (|1.9p and (|1.12p we have 

Sc{t)Qc - S{t)Q = f dudq ICc Sc{t)fc{p')Sc{t)fc{q') 

dujdq K-c Scit)fc{p)Scit)fciq) 
diodq /Coo S{t)f{p')S{t)f{q') 



+ / dojdq /Coo S{t)fip)S{t)f{q). (5.7) 



As usual, we split the integrals via 



1 - ^\p\ + \q\<r{c) + l|p| + |9|>r(c)- 



Here r(c) = (logc)"' with 1/2 < ai < 1. 

When the momentum are large we have the estimate 



R«xS"-i 



dujdqdp l|p| + |g|>r(c) l^c Sc{t)fc{p')Sc{t)fc{q') 

<bl j dwdqdp l|p|+|,|>^(c) J^{p')J^{q') 



<s«-i 



= bl j dwdqdp l|p|+|q|>r(c) K^c J'\p)j'^{q) 

< bljf"\r{c)) [ dudqdp l|p| + |,|>.(c) Jf"^{p)jf"\q) 

< Ablaojf"\r{c)) [ dudqdp l|p| + |,|>.(c) jf"\p)J^'H,q) 

JR«xS"-i 



< AJ^/2(r(c)) 



Above we have used the estimate for ICc{p, q,w) from Lemma 13.31 It is crucial that 
we gain back the square in the relativistic Maxwellian as in p.Sp , in order to use the 
estimate p. lip in the last line. We have also used the upper bound of J^[p) in the 
weight of the norm || • ||c from Theorem 12. 2| and the pre-post collisional invariance 
p.Sp . The independent of the speed of light bound for J{p) in Lemma [33] would not 
be enough. All the other terms on the r.h.s. of (|5.7p can be easily controlled by the 
same upper bound when the momentum are large. 

For bounded momentum we use more splitting. Notice that the momentum post- 
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collision depend upon the speed of light. We expand the terms on the r.h.s. of (|5.7p : 

dojdq l|p|+|g|<r(c) l^c Sc{t)fc{p)Sc{t)fc{q) 

dujdq l|p|+|g|<r(c) Sc{t)fc{p)Sc{t)fc{q) 
dujdq l|p|+|g|<r(c) ^oo S{t)f{p')S{t)f{q') 
dujdq l|p|+|g|<r(c) ^oo S{t)f{p)S{t)J{q) 



= I + II + III. 

For the first term we have 

dujdq l|p| + |,|<.(c) {ICc-lC^} {S{t)f{p')S{t)f{q')-S(t)f{p)S{t)f{q)}. 
For the second term in our splitting above is 



= j diodq l|p| + |q|<r(c) 


/Ce [S,{t)Up') - 


- Sc{t)f{p')]S,{t)Uq') 


(5.9) 


+ j dujdq l|p| + |q|<r(c) 


JCc [Sc{t)f{p') - 


S{t)f{p')]S,{t)Uq') 


(5.10) 


+ J dujdq l|p| + |g|<r(c) 


/Cc [S4t)Mq') - 


~Scit)f{q')]Sit)fip') 


(5.11) 


+ J dudq l|p| + |q|<r(c) 


ICc [Sc{t)f{q') - 


- Sit)fiq')]Sit)fip') 


(5.12) 


+ j dujdq l|p| + |g|<r(c 


) /C, [S{t)fiq') - 


- S{t)f{q')]Sit)fip') 


(5.13) 


+ J dujdq l|p| + |q|<r(c; 


, /Ce [Sit)fip') - 


S{t)f{p')]S{t)fiq'). 


(5.14) 



The third term above can be expanded as 

dujdq l|p|+|9|<r(c) ICc [Sc{t)fc{p) - Sc{t)f{p)]Sc{t)fc{q) 

dudq l|p|+|g|<r(c) [Sc{t)f{p) - S{t)f{p)]Sc{t)fc{q) 
dujdq l|p|+|5|<r(c) ICc [Sc{t)fc{q) - Sc{t)f{q)]S{t)f{p) 

dudq l|p|+|5|<r(c) [Sc{t)f{q) - S{t)f{q)]S{t)f{p). 



We will now estimate the first and second term in full detail, and we will then explain 
how the third term can be estimated in exactly the same way as the second term. 
To estimate term I, we use the following bound from Lemma 13.31 

\ICc-ICoo\<A{l + \p\ + \q\f/c\ 

For simplicity without loss of generality we have fixed m = n = 9. Thus we have 

Pllixroo < Ablricf/c' f dujdqdp /(p)/(g) < Ab'^ric)yc\ 
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This follows from the uniform estimate H/'^jloo !i ^2 assumed in Theorem [ 

Next we estimate each of the terms in II. With the Glassey-Strauss variables 
([THI)) . pTTTj) and when \p\ + \q\ < r(c) we have for some A>0 that 

b'l < b' - P'l + Ip'I < + 2 (H + M) . 

A similar estimate also holds for \q'\. We conclude that \p\ + \q\ < r(c) implies 

b'l + \q'\ < Aric). 

We will use this argument implicitly several times in the following estimates. 

For (|5.9|) . we use a pre-post coUisional change of variable and end up with an 
upper bound in terms of our unknown difference. After integration over dp, we have 

dudqdp l|p|+|g|<r(c) [Sc{t).fc{p') ~ Sc{t),f{p')]Sc{t)fc{q') 

< I dujdqdpl\p,\+\q>\<Ar(c) K^c WMp ) - f{p)\\L^\\Mq')\\L^ 

dojdqdp l|p|+|q|<Ar(c) ICc{p,q,(^) \\fc{p) - f{p)\\L^\\fc{q)\\L^ 

< Ac7oric)\\fc - /||liloc62. 



< 



As usual 62 is the constant from Theorem 12.21 We have also used the estimate p.2p 
combined with the exponential decay in J^{q) above. 

We look now at the difference (|5.10|) . which is a difference in the the spatial 
argument. We expand it out as follows 

S,{t)f{p') - S{t)f{p') = f{t,x- tp,p') - fit, X - tp,p') 

= f{t,x-tp + t{p-p),p) - f{t,x- tp, p) 

= T^S{t),f{p') - S{t)f{p'), u = t{p- p). 

On \p\ < r(c), with < < < T, and c > 1, we have 



\u\ = \t{p-p)\=t 



<L(^l + \p\^)<2^r{cf. (5.15) 



Taking the supremum in the x variable, integrating over dp and putting all of the 
estimates in this paragraph into (|5.10p . we obtain the upper bound for || (j5.10p || titcx, 

< J d^dqdp lH,|,<.(.) /C. \\SAt)fip') - Sit)fip')U^\\S.it)fM\Lr 

< sup / dujdqdp l\p,\+\g,\<Ar{c) ICc \\Tufip) - fip)\\L^\\fc{q')\\L^ 



< sup / rfwrfqdp l|p|+|,|<^r(c) /Cc - /(p)|1l-||/c('7)||l° 

\u\<j^ric)^J 

We have done a pre-post coUisional change of variable. Furthemore 

dp 



< A(7ob2r{c) sup / 

|„|< T r(c)3 J 



|«|<^r(c)3 

<Aaor{c) sup - /Ulilj 

|tt|<^r(c)3 
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These last estimates follow from the upper bound for in Theorem 12.21 combined 
with the estimate for JCc in p.2p . and also Lemma 

For (|5.11[) . we use the method from (|5.9p . Then for (|5.12[) . we use the method from 
(jS.lOp . except now we use the symmetric reverse estimate for the kernel /Cc(p, 9, w) as 
given in the display just below p.2p . 

For (|5.13|) . we estimate \q' ~ q'\ as in p.3|) . With that, we have 

S{t)f{q') - S{t)f{q') = S{t)f{q') - S{t)f{q' + [q' q')) 
^S{t)f{q')-TlS{t)f{q'), u = q'-q'. 

Now the rest of the estimate for (|5.13p can be handled exactly as in the estimate for 
(|5.10p , the only difference being that the difference is in terms of rather than . 
The end result is that 

\\mSI\\L^^L^<Aaor{c) sup \Kf - fh^L^. 

|«|<^r(c)3 

The remaining term (|5.14p needs to be estimated in a slightly different way. 

Lastly we estimate the term (|5.14p which forms a difference in the momentum 
arguments of the same function. For this term we expand the differences as 

S{t)f{p') - S{t)f{p') = S{t)f{p' + ip'-p')) - S{t)f{p') 
^TlS{t)f{p')^S{t)f{p'), u^p'-p'. 

On IpI + \q\ < r(c), with p.3|) . and c > 1, we have 

\u\^\p' -p'\<A 2 



r(c)3 



Taking the supremum in the x variable, integrating over dp and putting all of the 
estimates in this paragraph into (|5.14p , we obtain the upper bound 

liraiUiLsc < J dc^dgdp l|p| + |,|<,(,) /C, \\S{t)fip') - Sit)f{p')U^\\S{t)f{q')h^ 

< sup / dwdgrfp l|j,|+|,|<^(c) /Cc ||t?/(p') - /(p')IU~II/(^')IIl~ 



<<7Qr{cf/'^ sup / dujdqdpl\p\+\q\<r{c)\\Tlf{p) - f{p)\\L-^\\f{q)\\L-^ 
<ai)r{cf/'^ sup / dwdgdp l|p|+|q|<4^(c) |1t?/(p) - /(p)IIloo 

We have used (|3.2p . and the pre-post change in the Newtonian variables Then 
Wm^WLiL^ < Aaoricy/^ sup ||r?/ - /lU^^^^. 

This completes all of our estimates for the terms in II. 

Notice that each of the terms in III has an analogous term in II, except there are 
not post-collisional velocities to be concerned with. Thus the estimates in III follow 
exactly as the estimates in II with the same result. 
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We combine all of the estimates for each of the terms I, II, and III to obtain 

+ AMcV^' sup - fW^.L^ + IK/ - f\\^.^^) + a/-^. 

From this inequality and Lemma |5. II we conclude the statement of Lemma l5.2l for any 

small e > 0, since r(c) = (logc)"^ with 1/2 < ai < 1. □ 

Lastly, we prove the main Newtonian Limit Theorem 12.21 

Proof. We consider the difference of solutions to (|5.ip and (|5.2p with c > 1 as: 



fcit,x,p) - f{t,x,p) = Sc{t)fo,c{x,p) - S{t)foix,p) 

+ / ds Scit - s)Qcis,x,p) - S{t - s)Q{s,x,p). 
Jo 

After taking the L^L^ norm we obtain the following inequality 

+ [ ds\\S,{t-s)Q,{s)-S{t-s)Q{s)U.r^. 
Jo 

By translation invariance of the L'^ norm, we have 

\\Sc{t)fo.c - S{t)fo\\LlL^ < ||5,(t)/o,c - 5c(t)/o||LiLS= + \\Sc{t)fo - S(t)/o||LiLjo 

= ll/0,c - Ml^L^ + \KSit)fo - Sit) foil Ll,L^, 

with u = t{p — p). When \p\ < r(c), and r(c) = (logc)"^ with 1/2 < ai < 1, we use 
the estimate in (|5.15p to conclude 

\u\<AT/c'^-^, V£e(0,2]. 

Thus under the assumptions of Theorem 12.21 on \p\ < r{c), we have 

sup \KS{t)fo - S{t)M\L. rr < MT)/c'-', Ve G (0,2]. 

\u\<AT/c^-^ 

On the other hand, if \p\ > r(c), then we use estimates in the spirit of (|5.8p : 

\\Sc{t)fo - '5'(t)/o||Li(|p|>r(c))L;f < 2||/o||Li(|p|>r(c))L~ 

= 2 dp l|p|>r(c) II/o(p)||l~ < 2^2 / dp l|p|>r(c) (p) 

< A ^lP/\r{c)) f dp l|p|>,(,) 

< A^^^/\ric)) < A < ^. 



The estimates above follow from the bounds on the limit solution in Theorem 12.21 and 
Finally the non-linear asymptotic estimate from Lemma 15.21 vields 

\\SS-s)Q,is)-Sit-s)Q{s)\\LiL^ 

<Aric) \\f,{s)-f{s)Uir^+A{T)/c'-', We e [0,2]. 
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Putting all of these estimates together we have 

\\.m - fmLiLr < Whc - fohiL^ + A{T)ia- 

+ Ar(c) / |l/c(s)-/Gs)|Uiioo. 

^0 

By the Gronwall inequality 

ll/c(t) - /WIUiL5= < (||/o,c - hWLlL^ + A{T)I^-^) e^^'^(^). 

By the slow growth of r(c), with p.lOp . this gives easily the Cauchy property with 
the convergence rate of as stated in Theorem 12.21 □ 

This last proof completes the main focus of this paper. We conclude with two 
Appendicies. In Appendix A, we elucidate a collection of Lorentz Transformations 
that may be useful for choosing coordinates as in (|1.17p . In Appendix B, we point 
out some differential cross-sections which appear in the physics literature for the 
relativistic Boltzmann equation. 

Appendix A: Lorentz Transformations. In this appendix, we will write down 
three Lorentz Transformations which may be useful in choosing coordinates in rela- 
tivistic Kinetic Theory. In the prior sections, we worked in arbitrary dimensions 
N > 2. For simplicity, in this section we work in the physical dimension = 3. In 
this section we will also use the notation 

For more on Lorentz Transformations, we refer to |501l61j. We give their definition: 
Definition 5.3. A four by four matrix A is a Lorentz Transformation if 

(AP) • (AQ) = P-Q. 

Let D = diag{l, —1, —1, —1). Equivalently, A is a Lorentz Transformation if 

A^DA = D. 

Any Lorentz Transformation, A, is invertible and the last display implies 

A-i = DA^D. 



Thus the inverse given in (jl.l7l) can be easily computed. Since A has sixteen compo- 
nents and is restricted by ten equations, A has six free parameters. 

Perhaps the most well known Lorentz transformation is the boost matrix. Given 
V = {v^,v'^,v^) £ M.^, we write the boost matrix as 



Ah = 



7 



+ (7-1)^1^ 



(7-1)7^ 
l + (7-l)^ 



(7-1)^ 
l + (7-l)^i^ / 
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where 7 = (1 — Notice that A;, has only three free parameters. It is weh- 

known that any Lorcntz Transformation can be expressed as a rotation multiphed by 
a boost matrix. 

We are exclusively concerned with proper orthochronous Lorentz transformations 
A which send P + Q into the center-of-momentum system 



A(P + 0) = (V^, 0,0,0)* 



(5.16) 



We recall that s is defined in (|1.2p . We will give three examples of Lorentz Transfor- 
mations satisfying (|5.16p . 

Example 1: The Boost Matrix. We choose v such that A;, satisfies (|5.16p . Let 



p + q 

Po+qo' 



Then A;, is given by 



/ po+qo 
pi+gi 

P2+q2 



pq + go 



P2+22 



+ (7-1)^ 
(7-l)7¥ 



(7-1)7^ 
+ (7-1)t¥ 



(7- 



_ £3+23 

(7-in¥ 
(7-1)t¥ 



^ i + (7-i)^; 



where ^j^^ = ^ gy a direct calculation, this example satisfies (j5.16p . We 

notice that Ab — > /4, the four by four identity matrix, as c — > 00. 

Example 2. Given p,q E M.^ with p + q 0, we can always choose three orthonormal 
vectors 

...1 ...2 ,..3 P + q 



- ^ - ^ - \p + q\ 
If, for example, pi + qi ^ and P2 + q2 ^ 0, then we can explicitly write 

1 _ (-(P2 + q2){pz + 93), 2(pi + gi)(p3 + 93), ~{P2 + 92)(P3 + 93)) 



w 



I (-(P2 + g2)(p3 + 93), 2(pi + gi)(p3 + 93), -(P2 + g2)(p3 + 93)) I 

If instead pi + qi = 0, then 

= (1,0,0), w'^ = X /\w^ X w^\. 

With this notation, we can write down a Lorentz Transformation which maps p-\-q 
as 



A 



Po+qo 
\P+Q\ 





pi+qi 
\P+Q\ 
w\ 
wl 


P2+q2 
\P+Q\ 
W2 
wl 


P3+q3 
\P+Q\ 
wl 
wl 


\p+q\ 


pi+qi Po+qo 


P2+q2 Po+qo 


Pa+q-i Po+qo 


\P+Q\ 


\p+q\ \P+Q\ 


\p+q\ \P+Q\ 


\p+q\ \P+Q\ 



(5.17) 



This Lorentz transformation clearly satisfies ()5.16p . This example is fully derived in 
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Example 3: Hilbert- Schmidt form for the Linearized Operator. We now wish to derive 
a Lorcntz Transformation which satisfies (|5.16p but also 



A(P-O) = (0,0,0,5)* 



(5.18) 



In |151 p. 277] this transformation is used to write down a Hilbert-Schmidt form for the 
hnearized collision operator. But the transformation is not written down explicitly 
therein. What is given is conditions (|5.16p . (|5.18p and the first row Aq. From this 
information we can make the following educated guess for A satisfying (|5.16|) and 
(EUl): 



/ Po+go 



A: 



VP— go 

a 



pi+qi 


P2+92 


P3+g3 








A} 


A^ 


AJ 


(PX(})1 


(pX<j)2 


(pxg)3 


|pxg| 


|pxg| 


Ipxrjl 




P2 — 92 


P3 — 93 


a 


g 


a 



By the Lorcntz condition in Definition 15. 3[ we may conclude 



AJ 



\p X q\ 



2bx g| ^ 

^/{P~Q? 



The coefficients A^ [i = 1,2,3) are similarly completely determined by the Lorentz 
condition. However, this turns out to be a long computation. In the end we obtain 



^ 



2 {p^ {pa - qpP ■Q} + q^ {go - PqP ■ Q}) 
VsglP X q\ 



1,2,3). 



Now, we have a complete description of this Lorentz transformation in terms of p, q. 

The way we derived this matrix was by guessing the first, third and fourth rows 
and then using the Lorentz condition. However, it is much easier to verify that this is 
indeed a Lorentz Transformation after the fact by checking the condition in Definition 
15.31 The interested reader may find this calculation in [55] . 

Appendix B: Examples of relativistic Boltzmann cross-sections. In [?Tj . 

Grad gave several examples of important collision kernels to study in the classical 
theory. The same type of discussion would be important for the relativistic kinetic 
theory, in this direction we point out [19] . Here we write down a few scattering kernels 
which can be found in the physics literature. 

The calculation of the differential cross-section in the relativistic situation utilzes 
quantum field theory, see for instance |51j . Cross sections are not presently derived 
from a scattering problem (as it is done in the Newtonian case) in particular because 
there is no widely accepted theory of relativistic N-Body dynamics at present. 

Above, we have written everything down with the mass m normalized to one, 
m = 1. Here we write down the mass without normalization. 



Short Range Interactions. [25j[5T]. For short range interactions, 



(T = constant. 



(5.19) 



This "hard-ball" cross section is the relativistic analouge of the hard-sphere kernel in 
the Newtonian ease. Indeed, as we have already mentioned, the Newtonian limit of the 
relativistic Boltzmann equation in this case is the hard-sphere Boltzmann equation. 
This is also sometimes written as scr — constant. 
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M0ller Scattering. [131 P-350]. M0ller scattering is used as an approximation 
of electron-electron scatttcring. 



where the magnitude of total four-momentum scaled w.r.t. total mass is 

2mc' 

2 

and ro = ^^^^a is the classical electron radius. Photon-photon scattering is often 
neglected because the size of the cross-section is 'negligible'. 

Compton Scattering. [15l p. 351]. Compton scattering is an approximation of 
photon-electron scattering. 

1.2/1 , 1 ^^(l-cos6l)2 p-(l-iO(l-cos( 



2"' ^'1^ 41-i^(l-cos6i) V l-5^(l-cos( 
where 

2 2 

^ = 1 . 



See [251 P-81] for the Newtonian limit in this case. 

(elastic) Neutrino Gas. [I7l p. 478] For a neutrino gas, the differential cross 
section is independent of the scattering a angle 0. 

q2 

where G is the weak coupling constant and 27rfi, is Plank's constant. These are massless 
particles. See also [TSl p. 290]. 

Israel particles. [40l p. 1173] The Isreal particles are the analogue of the well 
known "maxwell molecules" cross section in the Newtonian theory. They are 

m b{0) 

^[9,0) - 



2.9 1 + {g/mcY 

With this cross section, Israel derived eigenfunctions for the Linearized relativistic 
Boltzmann collision operator. Variants are used in [T31[S2]- For instance the cross 
section for "Maxwell Particles" is formed by removing the factor 1 + [g/mcf . It 
converges to the maxwell molecules cross section in the Newtonian limit. 
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